Why Dirac points in graphene are where they are? 
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We present a simple group theory explanation of the fact that the energy bands merge in the 
corners of the Brillouin zone for graphene and for two particular cases of Kagome lattice for arbitrary 
tight-binding Hamiltonian. We connect the linearity of the spectrum in the vicinity of these points 
for monolayer graphene, bilayer graphene for AA stacking and Kagome lattice with the properties 
of the conical points of the surface, known from geometry. 



INTRODUCTION 

The contacts of energy bands have been studied from 
the early days of quantum sohd state physics P, Q • The 
early classical papers in the field are reprinted in Ref. Q . 
The existence of such points of contact in graphene was 
also known since the early quantum mechanical studies 
of graphite. When graphene was first isolated experi- 
mentally [3] , these points were immediately brought into 
the focus of attention of both theorists and experimental- 
ists. The reason is very simple: for undopcd and ungated 
graphene these points lie on the Fermi surface, and many 
physical properties of such graphene are determined by 
the electron states in the vicinity of these points. In addi- 
tion, the electron spectrum in the vicinity of these points 
is linear (in a monolayer graphene), and these states are 
effectively described by the Dirac equation [5i]. This is 
why these points of contact in monolayer graphene (and 
other materials showing similar properties) are called the 
Dirac points. 

Graphene is traditionally considered within the tight- 
binding approximation. The approximation starts from 
mentioning that in isolated form, carbon has six electrons 
in the orbital configuration \s^2s^2p^ . When arranged 
in the honeycomb crystal, two electrons remain in the 
core Is orbital, while the other orbitals hybridize, form- 
ing three sp^ bonds and one pz orbital. The sp^ orbitals 
form the a band, which contains three localized electrons. 
The bonding configuration among the Pz orbitals of dif- 
ferent lattice sites generates a valence band, or 7r-band, 
containing one electron, whereas the antibonding con- 
figuration generates the conduction band (tt*), which is 
empty 0. 

The early classical papers consider merging of Bril- 
louin zones 0. On the other hand, in the framework of 
the tight-binding approximation there exists only a single 
Brillouin zone. However, this zone contains the number 
of energy bands, equal to the number of atoms in the 
elementary cell (more exactly, the number of orbitals per 
elementary cell considered), which emulate a finite num- 
ber of Brillouin zones. 

The structure of graphene can be seen as a triangu- 
lar lattice with a basis of two atoms per unit cell. For 
a model tight-binding Hamiltonian, taking into account 



only nearest-neighbor and next-nearest-neighbor hop- 
ping, it can be easily shown that the energy bands these 
bands neither overlap, nor have a finite gap between 
them, but touch each other in two points in momentum 
space, both for monolayer and bilayer graphene. We 
show in the framework of group theory that taking into 
account hopping between more distant neighbors does 
not shift these points of contact; their positions are de- 
termined only by the point symmetry of the honeycomb 
lattice. The linearity of the spectrum of the electrons 
in the vicinity of the points of contact in monolayer and 
bilayer graphene with the AA stacking also remains for 
the general tight-binding Hamiltonian. 

In addition we prove the existence of the Dirac points 
in the corners of the bands for the isotropic Kagome and 
Lieb lattices, described in the framework of the general 
tight-binding Hamiltonian. 

THE GENERAL TIGHT-BINDING 
HAMILTONIAN 

The structure of graphene can be seen as a triangular 
lattice with a basis of two atoms per unit cell, displaced 
from each other by any one (fixed) vector connecting two 
sites of different sub-lattices, say 5 = —a (1, 0). The gen- 
eral Hamiltonian for the tt bands is 

VEa^*(a + '5)e-*(a+^) S^t'(a)e'k.a J, 

(1) 

where a is an arbitrary lattice vector, that is a linear 
combination of ai = | (3, \/3) , a2 = f (3, — \/3) . 
The selection rule for matrix elements [6j] gives 

^t(a + ,5)e''^-('^+'') =0, (2) 

a 

where K is a the corner of the Brillouin zone. In fact, we 
are dealing with the product of two functions: i(a -I- 5) 
realizes the unit representation of the point symmetry 
group C3 (the full symmetry group of the inter-sublattice 
hopping is Ca^, but the restricted symmetry C3 will be 
enough to prove the cancelation) . As far as the function 
giK (a+(5) concerned, rotation of the lattice by the an- 
gle 27r/3, say anticlockwise, is equivalent to rotation of 



2 



the vector K in the opposite direction, that is to sub- 
stitution of the three equivalent corners of the Brilluoin 
zone: Ki ^ Ka ^ K3 ^ Ki, where Ki = (|f , 



K2 - 



and Ka 



_ 27r 
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to the rotation g*^' is muhiphed by the factor 
and reahzes x — iy representation of the group C3. Be- 
cause each of multiphers in Eq. ^ reahzes different ir- 
reducible representation of the symmetry group, the ma- 
trix element is equal to zero. Simply speaking, at a point 
K the sublattices become decoupled, and this explains 
the degeneracy of the electron states in this point (these 
points) or, in other words, merging of the two branches 
of the single Brilouin zone. 

On the other hand, generally 



Thus due 



(3) 



To understand this statement consider the maximum 
symmetry group of the intra-lattice hopping: C'ey. The 
function t'{a.) realizes the Ai representation of the group. 
Applying Eq. (??) we see that reducible representation 
of the group C^v , realized by the two functions e'^ ** and 
e'^ can be decomposed as Ai + B2. 

In addition, the tight-binding model provides us with a 
simple explanation why the dispersion law in the vicinity 
of the merging points is linear, that is why these points 
are Dirac points. The dispersion law for the Hamiltonian 
([T]) is given by equation 



F(£;,k) = 0, 



(4) 



where 



F{E,k) 



E + Ea i'(a)e*-=^ t(a + ,5)e'k.(a+5) 

t*{a + 5)e-''' ('^+*') E + J2j'i^)e^" 



(5) 



a 



dF 

dE 



0, 



dF 
5k 



0. 



(6) 



Recalling the rule for differentiating of a determinant, 
we realize that Eq. ^ guaranties that the conditions © 
for k = K. This explains linearity of the spectrum in the 
vicinity of the points K(K'). The axis of the cone is per- 
pendicular to the kx, ky plane, and the cone is a circular 
one, because any vector in the kx,ky plane compatible 
with the symmetry C31, is identically equal to zero, and 
any tensor of rank two compatible with the symmetry is 
proportional to the unity tensor. 

The important role played by discrete symmetries in 
protecting a fc-linear dispersion in graphene was pointed 
out in Ref. 0]. The appearance of massless Dirac 
fermions under conditions of hexagonal symmetry was 
considered in Ref. ■ Group theory was used to derive 
an invariant expansion of the Hamiltonian for electron 
states near the K points of the graphene Brillouin zone 
in Ref. [ll|. 



BILAYER GRAPHENE 

The same selection rule for matrix elements, we applied 
for the case of monolayer graphene, can be applied to 
bilayer graphene. 




Figure 1: Top view of a graphene bilayer; white and 
black circles: top layer carbon atoms; gray and red: 
bottom layer. [Adapted from Ref. [13j ] 



In mathematics the Dirac points, we are dealing with, 
are called conical points of the surface; if the surface is 
given by Eq. (|3|), the conditions for the conical points stacking has the form 



The general tight-binding Hamiltonian for the AB 



H 



„«ka 



Eai4(a + <5)e«k.(a+5) 
Ea*'(a)e*'^ 



Eai3(a 
Eai4(a 
Ea^(a- 

Ea^'C 



. J)g«k.(a+5) 



a e 



(7) 



where t' and t are the intra- and inter-sublattice intra- 



layer hopping integrals respectively, and ti, and are 



3 



the A— A, B — B, and A—B inter-layer hopping integrals zero, due to symmetry reasons mentioned above. Hence, 
respectively at the corner of the Brillouin zone the Hamiltonian ([7]) 

For k = K(K)', the terms containing t,t^ and t4 are takes the form 



H 



/ Ea*'(a)e''^" 

\ Eai'(a)e'^" 



(8) 



We again have merging of the energy band. Of course, in The general tight-binding Hamiltonian for the AA 
this case the points of contact are not the Dirac points. stacking has the form 

I 



Eai*(a + '5)e-*'^-(=^+*) Ea^'(a)e''"" Ea ^4(a + <5)e''"('^+^) Ea ^i(a)e^'"" 

Ea^i(a)e*'^-=^ Eai|(a + <^)e-*-'"+^' Ea ^'(a)e^'"" Zji^ + S)e'''<''+^^ 



where ti and are the A — A{B — B) and A — B inter- layer hopping integrals respectively. At the corner of the 

Brillouin zone the Hamiltonian © takes the form 



/Ea^'(a)e'^-^ Eaii(a)e'^-^ \ 

Ea*'(a)e''^" Ea^i(a)e''^" 

Eaii(a)e'^- Ea^'(a)e'^'^ 

V Ea^i(a)e^^=^ Eai'(a)e'^V 

(10) 



and we again have merging of two bands. However, for 
the AA stacking, in distinction from the AB stacking the 
function 

F{E,\^) = \EI-H\, (11) 

where / is the unity matrix, satisfies Eq. Hence the 
points of contact are the Dirac points. 



KAGOME LATTICE 

Kagomc lattice has three sublattices. The general 
anisotropic tight-binding Hamiltonian for such lattice is 
characterized by three sets of exchange integrals. We con- 
sider general tight-binding Hamiltonian for two particu- 
lar cases. First consider isotropic case ti — t2 = <3 = t. 
The Hamiltonian in this case is given by equation 



H ^ 



Eji^y^"" X;a*(a + ai/2)e''^-('^+=^i/2) J2j{a + a3/2)e'^<''+''^^^^ 

X;^^* (a + ai/2)e-'i'-('*+'^i/2) Ea ^(a)^'''''' i (a + a2/2) e^'^-^^+^^Z^) 

t* (a + a3/2) e-»k.(a+a3/2) ^* ^ jj^/2) e-»k.(a+a./2) t(a)e'''"^ 



(12) 



where ai and a2 are the elementary lattice vectors, and Eq. (fT^ can be presented as 

^'=^'-^'- I y(k) r(k) eos (iifi) r(k) cos {^) \ 

r(k)cos(^^) y(k) T(k)cos(^^) , 

\r(k)cos(^^) T(k)cos(i^) Vi\.) I 
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(13) 



with obvious expressions for y(k) and T(k). We see, that 
the dispersion law for a general tight -binding Haniilto- 
nian can be presented as 



E{k) ^V(k)+T{k)E{k) 



(14) 



where E{]i.) is the dispersion law given by the Hamilto- 
nian taking into account only the nearest-neighbor hop- 
ping [13 



H = 




(15) 



One branch of the spectrum we find by inspection: 

^i(k) = -l. (16) 

The other two branches are given by the reduced Hamil- 
tonian 



1 



y3 



1 



(17) 



At the point K(K') the non-diagonal matrix elements of 
the Hamiltonian (fTT]) are equal to zero due to the reasons 
presented above, and we have merging of the bands. In 
addition, at the point K(K') the function F, defined by 



F{E,k) 



l~E 



l-E 



satisfies conditions ([6|). Thus we have two Dirac points 
in the corners of the Brillouin zone. Note also, that at T 
point one of the branches given by Eq. ([T7]) touches the 
branch Ei. 

Now consider particular case of the Kagome lattice 
ti = t2 — t, — 0. In this case, called isotropic Lieb 
lattice, we recover Eq. ([T^ , the role of the Hamiltonian 
(fT5)l being played by the Hamiltonian 







H 





cos (^) 



cos 







^ . (19) 



The spectrum is 



i^i(k) = 0, i^|_3(k) = cos^ 



k • ai 



cos 



k • ai 



(20) 



This time all three bands merge at the point M = 
(|j,0), the zones E2 and E3 again forming Dirac cones 
in the vicinity of the point of contact. 

In this work we presented simple arguments explain- 
ing merging of the energy bands in the corners of the 
Brillouin zone for the case of both monolayer and bilayer 
graphene and two particular cases of the Kagome lattice 
and the Dirac nature of these merging point in the case 
of monolayer graphene, bilayer graphene for AA stacking 
and Kagome lattice. 
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